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TECHNICAL PAPER 


FURTHER DEVELOPMENTS IN MODELING DIGITAL CONTROL SYSTEMS 
WITH MA-PREFILTERED MEASUREMENTS 

I. INTRODUCTION 


Books on modern digital control systems usually treat the problem of controlling a linear time- 
invariant plant, driven by a zero-order-hold with a sampled output, as shown in Figure 1 (see, for 
example, Reference 1). The state vector x(t) € R n ; the control input vector u(kT) € R'; the output or 
measurement vector yi(kT) 6 R m ; the system matrix F € R nxn ; the control matrix G € R nxr ; the output 
matrix Cj € R mxn . It is well known that this system can be modeled at the sampling instants by the discrete 
state equations [1] 


x(k + IT) = A x(kT) + B u(kT) 

(1) 

y,(kT) = C, x(kT) 

(2) 


where 


4>(t) =*zf' [(si-F)- 1 ] e R nxn , 

(3) 

A = 4>(T) € R nxn , 

(4) 

B = [/ 4>(\)d\] G € R nxr . 

(5) 


o 


In light of equation (2), yi(kT) represents an instantaneous measure of the system at the sampling instant 
kT and, hence, can be thought of as being an instantaneous measurement vector. 

Unfortunately, not all systems in the real world produce discrete measurements like this. For 
example, there exist systems where discrete measurements represent average measures of the system 
over the time interval T between outputs. Such systems can be found in the aerospace field wherever 
startrackers and some state-of-the-art rate gyroscopes annd accelerometers are used, to name a few [2]. In 
one system of this type, discrete measurements are generated every T/N seconds and averaged as in 
Figure 2. Every N measurements are averaged to produce the averaged measurement vector y A (kT) € R h , 
every T seconds. The system in Figure 2 also allows for the possibility of instantaneous measurements. 
Discrete state variable models for this system were derived in References 2 and 3. 















The system in Figure 2 was then generalized to that shown in Figure 3 and discrete state variable 
models for it were derived in Reference 4. In Figure 3, the output matrix C F € R pxn and so the output 
vector z(t) 6 R p . This is sampled every T/N seconds and multiplied by the weighting matrices Hj € R qxp , 
j = 0,1 ,...,N-1 . After N repetitions, the results are summed to produce the weighted-averaged mea- 
surement vector y F (kT) € R q every T seconds. Functionally, this is equivalent to passing the mea- 
surements generated every T/N seconds through a multi-input/multi-output moving-average (MA) 
process with coefficient matrices Hj, j = 0, 1 1 [5]. The output of the MA prefilter is sampled every 

T seconds to produce the output vector y F (kT). Hence, y F (kT) can also be thought of as an MA-prefiltered 
measurement vector. Observe that the system in Figure 2 is a special case of the one in Figure 3. 

Reference 4 also derived discrete state equations for the system in Figure 4 when the matrix E_ € 
R qxr is uniquely chosen. The system in Figure 4 is the same as the one in Figure 3 except the MA- 
prefiltered measurement vector y F (kT) is modified to yield the output vector y F '(kT) € R q . This system is 
the basis for a new type state reconstructor that is developed in References 6 to 10. 

Enlightened by the work of Hagiwara and Araki [11], it became clear that the results for the sys- 
tems in Figures 3 and 4 could be extended to the systems shown in Figures 5 and 6. In Figure 5, the 
sampling intervals on the measurements to be MA-prefiltered vary from T/2 to T/N with the length of the 
MA prefilters varying from 2 to N. For i = 2,3,...,N, the output matrix C F j € R p,xn ; the output vector z ; (t) 
€ R Pl ; the weighting matrix Hj j € R q,xpi where j = 0, 1 ,.. ,,N-1 ; the MA-prefiltered measurement vector 
y Fi (kT) € R q '. Figure 6 is like Figure 5 except that each MA-prefilter measurement vector y Fi (kT) is 
modified to produce the output vector y F j'(kT) € R q '. The matrix Ej_ € R q,xr . Observe that Figures 3 and 4 
are special cases of Figures 5 and 6, respectively. 

This paper presents state variable representations for the systems in Figures 5 and 6. This is done 
in Section III by extending the results for Figures 3 and 4 that are summarized in Section II. Concluding 
remarks are made in Section IV. 


II. PREVIOUS STATE VARIABLE REPRESENTATIONS FOR DIGITAL 
CONTROL SYSTEMS WITH MA-PREFILTERED MEASUREMENTS 


In Reference 4, two discrete state variable representations are derived for the system in Figure 3. 
The first one is 


x(k+ IT) 


'a o" 


"x(kT)~ 
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p(k+ IT) 


D + 0 


p(kT) 


_E + 


u(kT) 


( 6 ) 
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C, 
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’ x(kT)" 

y F (kT) 
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_t](kT)_ 


( 7 ) 
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Figure 6. The system in Figure 5 














where the matrices A and B are defined by equations (3) to (5), the vector -q(kT) € R q , and the matrix I is 
a qxq identity matrix. The matrix D+ is defined by the relationship 


D+ = Ha + € R qxn 


where 


H = [H 0 j H, ]... ! H N .|] € R qx(Np) 


and 


C F <(>(t) 

C F <J>[T-(T/N)] 


a+ = 


£ J^(Np)xn 


C F <f>[T-(N-l)T/N] 


The matrix E + is defined by the relationship 


E + = H (3+ € R qxr 


where 


P + = 


T 

C F [f 4>(X)d\] G 


T-T/N 


C F [f 4>(X)d\] G 


T-(N-I)T/N 

C F [J 4)(\)dX] G 


£ ^(Np)xr 


The second state variable representation for the system in Figure 3, derived in Reference 4, is 
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u(kT) 


( 8 ) 


x(k+ IT) 


A 

o" 


x(kT) 


"b" 


= 
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T|(k+ IT) 


0 

0 


H(kT)_ 


I 


"yi(kT)“ 


“c, 

0 ' 


'x(kT)' 

y F (kT) 


D_ 

E_ 


H(kT) 


(9) 


where r|(kT) € R r and the matrix I is an rxr identity matrix. The matrix D_ is defined by the relationship 


= Ha_ € R qxn 


where 


a_ — 


C F <K0) 

C F <J> (-T/N) 


LCf 4) [-(N-I)T/N]J 


£ j^(Np)xn 


The matrix E_ is defined by the relationship 


E_ = H(3_ € R qxr 


where 


0 - = 


C F l/ 4>(X)dX] G 


r -(T/N) 

C F [J 4*(X)dX] G 
o 


-(N-IXT/N) 

C F [/ <KMd\] G 
o 


£ j^(Np)xr 


I 
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In Reference 4, the system in Figure 3 was modified to become that in Figure 4. A discrete state 
variable representation for it was shown to be 


x(k+ IT) = Ax(kT) + Bu(kT) 


y.(kT) C, 

yp'(kT) D 


x(kT) 


The system in Figure 4 and this discrete state variable representation for it are the basis for the new type 
state reconstructor developed in References 6 to 10. 

III. NEW STATE VARIABLE REPRESENTATIONS FOR DIGITAL CONTROL 
SYSTEMS WITH MA-PREFILTERED MEASUREMENTS 


Extending the results in Section II for the system in Figure 3, two discrete state variable represen- 
tations can be defined for the system in Figure 5. From equations (6) and (7), the first one is 


x(k + IT) A 

r| 2 (F+lT) d 2 . 


T] N (k + IT) Dp, 



x(kT) 

y\i(kT) 


T|N(kT) 


u(kT) 


yi(kT) 

YF2(kT) 


y FN (kT) 


nxH 

[•]N> 


x(kT) 

^l 2 (kT) 


In T}N(kT) 


where i = 2,3,- -.,N, the vector T|i(kT) € R q ', and the matrix I; isaqjXqj identity matrix. The matrix D i+ is 
defined by the relationship 
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Dj + = H ( a i+ € R q,xn 


where 


[H i0 ! H n | 


• • ! H i(i _ )) ] € R qiX<ip,) 


Cpi <J>(T) 


C Fi c}>(T - -L) 

1 


£ R(ip,)xn 


C Fi 4>[T - (i-1) -y-] 


The matrix E i+ is defined by the relationship 


Ej + = Hi3i+ € R 4 '’ 


where 


C Fi [/ 4»(X)dX] G 


C Fi [/ <j>(\)d\] G 


£ R(ip.)xr 


T-(i-l )T/i 

Cpj [J <f>(A)d\] G 

o 


From equations (8) and (9), the second state variable representation for the system in Figure 5 is 


x(k+ IT) 
T)(k+ IT) 


‘ [A 01[ x(kT)] p“ 

= + u(kT) 

J Lo o J L ri(kT) J Li. 
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y.(kT) 


c, 

0 

yF2(kT) 

= 

d 2 _ 

e 2 _ 

y FN (kT) 
Lr _j 


D n _ 

En- 


- x(kT)‘ 
_ 3(kT) 


(15) 


where all vectors and matrices have been previously defined except for the following. For i = 2,...,N, the 
matrix Dj_ is defined by the relationship 


Dj_ - Hi«j_ € R q,xn 


where 


Cpj d>(0) 


«i_ 


Cpj <t> ( — T/i) 


£ R(ip.)xn 


C Fi <J>Hi-l) T/i] 


The matrix Ej_ is defined by the relationship 
Ej_ = Hi(3i_ € R q,xr 


where 
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Pi- = 


C Fi [/ <J)(\)d\] G 


-T/i 


C Fi [f <J>(A)dA] G 


- — (i— l KT/i) 

^Fi [J <f>(A)dA] G 

o 


£ R(ip,)xr 


For the system in Figure 6, a discrete state variable representation can be formed by extending the 
results in Section II for the system in Figure 4. From equations (10) and (11), the result is 


x(k + IT) = Ax(kT) + Bu(kT) 


(16) 


y.(kT) ~ 


"c, " 

YF2'(kT) 

= 

d 2 - 

y FN '(kT) 


Dn- 


x(kT) 


(17) 


where all vectors and matrices have been previously defined. 


IV. CONCLUDING REMARKS 


This paper presented new general state variable representations for digital control systems with 
MA-prefiltered measurements. The previous models, derived in Reference 4 for the systems in Figures 3 
and 4, were extended to the systems in Figures 5 and 6. 

For the system in Figure 5, two models were presented. One is defined by equations (12) and (13) 
while the other by equations (14) and (15). As to which one is best to use for a given system, the follow- 
ing advice is offered. Choose the one which yields the least number of states. The resulting discrete state 
equations are more likely to be observable and controllable. 

For the system in Figure 6, the model presented is defined by equations (16) and (17). It is 
especially intriguing because it has the same plant equation as the state variable model given by equations 
(1) and (2) for the system in Figure 1. 
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